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ABSTRACT 

We present a systematic spectral analysis with Suzaku of six AGNs detected in the Swift/BAT hard X-ray 
(15-200 keV) survey, Swift J0138.6-4001, J0255.2-001 1, J0350.1-5019, J0505.7-2348, J0601.9-8636, and 
J1628. 1-5145. This is considered to be a representative sample of new AGNs without X-ray spectral informa- 
tion before the BAT survey. We find that the 0.5-200 keV spectra of these sources can be uniformly fit with 
a base model consisting of heavily absorbed (log A^h >23.5 cm" 2 ) transmitted components, scattered lights, a 
reflection component, and an iron- K emission line. T here are two distinct groups, three "new type" AGNs (in- 
cluding the two sources reported by lUeda et al.ll2007l) with an extremely small scattered fraction (/ sca t < 0.5%) 
and strong reflection component (R = Vl/2tt > 0.8 where Q, is the solid angle of the reflector), and three "clas- 
sical type" ones with / scat > 0.5% and R < 0.8. The spectral parameters suggest that the new type has an 
optically thick torus for Thomson scattering (A^h ~ 10 25 cm" 2 ) with a small opening angle 9 ~ 20° viewed in a 
rather face-on geometry, while the classical type has a thin torus (A^h ~ 10 23 " 24 cm" 2 ) with 9 > 30°. We infer 
that a significant number of new type AGNs with an edge-on view is missing in the current all-sky hard X-ray 
surveys. 

Subject headings: galaxies: active — gamma rays: observations — X-rays: galaxies — X-rays: general 



1. INTRODUCTION 

The growth history of supermassive black holes (SMBHs) 
in galactic centers is a key question for understanding cosmic 
history. This can be traced by observations of Active Galactic 
Nuclei (AGNs), where mass accretion onto SMBHs is taking 
place. Many observations indicate that obscured AGNs are 
the major AGN population both in the local universe an d at 
intermediate-to-high redshifts (see e.g., Hasinger 2008). In 
fact, much of the SMBH growth is theoretically expected to 
happe n during an obscure d phase in the evolution of AGNs 
(e.g., [Hopkins et al. 2006). Thus, detecting all the obscured 
AGNs in the universe and revealing their nature is a funda- 
mental issue to be addressed by modern astrophysics. Our 
knowledge of this population is still quite limited even in the 
local universe, due the difficulty of detecting highly absorbed 
objects in most wavelength bands. 

Obscured AGNs, often referred as type 2 AGNs, are char- 
acterized by heavily absorbed X-ray spectra and/or by the ab- 
sence of broad emission lines in the optical spectra. The uni- 
fied model of AGNs attributes these differences to the block- 
age of the line of sight by a dusty torus. Recent X-ray and IR 
observations have revealed populations of AG Ns without any 
signa ture of AGNs in the optical band (e.g.. iMaiolino et all 
2003). This fact implies that AGN surveys relying on optical 
emission signatures (such as the [O III] A5007 line) may be 
incomplete. 

Sensitive hard X-ray observations above 10 keV provide 
us with the least biased AGN samples in the local universe, 
due to the strong penetrating power against photo-electric ab- 
sorption, except for heavily Compton thick (column density 
of log A^h > 24.5 cm" 2 ) objects. R ecent all sky hard X -ray 
surveys performed with Swift/BAT dTueller et alJ l2008h and 
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INTEGRAL dBassani et alj|2006t iKrivonos etaT]|2007h have 
started to detect many local AGNs with about 10 times better 
sensitivity than that of previous missions in this energy band. 
They have shown that roughly half of all local AGN are in- 
deed obscured AGNs. These hard X-ray samples contain a 
number of new AGNs, for which detailed follow-up studies 
have not been done yet. Some of these sources objects were 
not recognized to be AGNs previously even though there were 
optical imaging and spectroscopy data available. 

Our team has started to make follow-up observations of se- 
lected Swift/BAT AGNs with the Suzaku satellite to obtain the 
best quality broad band high energy data. The major purposes 
are (1) to discover new populations of AGNs, (2) to measure 
their spectral properties, which make it possible to construct 
the A^h distribution of a complete Swift/BAT sample, and (3) to 
accurately constrain the strength of the reflection component, 
which is an important yet uncertain parameter in population 
synthesis models of the X-ray background in exp laining its 
broad-band shape (lUeda et alJl2003b iGilli et al.ll2007l) . In par- 
ticular, the simultaneous coverage over the 0.2-70 keV band 
with Suzaku is crucial to discriminate degeneracy of spectral 
models under possible time variability and to determine the 
amounts of scattered and reflected components, key informa- 
tion to constrain the geometry of surrounding matter around 
the nucleus. 

The first results obtained from the Suzaku follow-up of two 
Swift/BAT AGN s led to the disc overy of a new type of deeply 
buried AGNs dUeda et alj|2007l) . Here we present the results 
of a detailed analysis of six Swift/BAT AGNs, including the 
two sources (Swift J0138.6-4001 and Swift J0601.9-8636) 
already reported there 4 . This paper constitutes the first one in 
the series of papers to be published on Suzaku-Swift AGNs, 
followed by the second paper by Winter et al. (in prepara- 

4 For these two sources, the results of the present paper, which are based 
on more recent calibrati on of t he ins truments and background estimation, 
supersedes those by Ueda et al. 1 2007), although the essence is not changed. 
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TABLE 1 
List of Targets 



SWIFT Optical/IR Identification R.A. (J2000) Dec. (J2000) Redshift Classification 



J0138.6^K)01 


ESO 297-G018 


01 38 37.16 


-40 00 41.1 


0.0252 


Seyfert 2 


J0255.2-0011 


NGC 1 142 


02 55 12.196 


-00 11 0.81 


0.0288 


Seyfert 2 


J0350.1-5019 


2MASX J03502377-5018354 


03 50 23.77 


-50 18 35.7 


0.036 


Galaxy 


J0505.7-2348 


2MASX J05054575-235 1 139 


05 05 45.73 


-23 51 14.0 


0.0350 


Seyfert 2 


J0601. 9-8636 


ESO 005-G004 


06 05 41.63 


-86 37 54.7 


0.0062 


Galaxy 


J1628.1+5145 


Mrk 1498 


16 28 4.06 


+51 46 31.4 


0.0547 


Seyfert 1.9 



NOTE. — The position, redshift, and classification for each source is taken from the NASA/IPAC Extra- 
galactic Database. 



tion). § 2 summarizes the sample and observations. In § 3, 
we mainly present the results of detailed spectral analysis with 
Suzaku. Discussion is given in § 4. We adopt the cosmolog- 
ical parameters (H , O m , Q\) = (70 km s" 1 Mpc" 1 , 0.3, 0.7) 
throughout the paper. 

2. OBSERVATIONS 

We observed six SwiftlBKY AGNs with Suzaku between 
2006 April and 2007 January in the AO-1 phase. The tar- 
gets are Swift J0138.6-4001, J0255.2-0011, J0350.1-5019, 
J0505.7-2348, J0601.9-8636, and J1628. 1-5145, whose 
identification in the optical or infrared band is ESO 297- 
G018, NGC 1142, 2MASX J03502377-5018354, 2MASX 
J0505 4575-235 1 1 39, ESQ 005-G004, and Mrk 1498, respec- 
tively dTueller et alj [2008). The basic information of our tar- 
gets is summarized in Table Q] 2MASX J03502377-5018354 
and ESO 005-G004 had not been identified as active galax- 
ies in the optical band, and revealed to contain an AGN for 
the first time with the detection of hard X-rays by the BAT. 
The sources were selected as the Suzaku targets because at 
the time of our proposals these "new" Swift AGNs had no 
spectral information available below 10 keV. Recently short 
follow-up observations of these targets have been carried out 
with XMM-Newton or the X-Ray Telescope (XRT) onboard 
Swift (IWinter et al.ll2008al) . 

Suz aku, the fifth Japanese X-ray satellite dMitsuda et alj 
l2007h . carries four X-ray CCD cameras called the X-ray 
Imaging Spectrome ter (XIS-0, XIS-1, XIS-2, and XIS-3; 
iKovama et al. 2007) as focal plane imager of four X-ray tele- 
scopes, and a no n-imaging instrument called the Hard X-ray 
Detector (HXD; iTakahashi etaT1l2007h consisting of Si PIN 
photo-diodes and GSO scintillation counters. XIS-0, XIS-2, 
and XIS-3 are front-side illuminated CCDs (FI-XISs), while 
XIS-1 is the back-side illuminated one (BI-XIS). 

In this paper, we analyze the data of the XISs and the 
HXD/PIN (HXD nominal position), which covers the energy 
band of 0.2-12 keV and 10-60 keV, respectively, since the 
fluxes of our targets above 50 keV are too faint to be detected 
with the HXD/GSO. Table Q] shows the log of the observa- 
tions. The net exposure of Swift J0255. 2-0011 and Swift 
J0505. 7-2348 is about 100 ks and 80 ks, respectively, while 
that for the rest 4 targets is 20 ks each 5 . Because XIS-2 be - 
came unoperatable on 2007 November 7 dDotani et al.ll2007l) . 
no XIS-2 data are available for Swift J0255.2-001 1 and Swift 
J0350. 1-5019. For the observations of Swift J0255.2-0011 

5 The primary goal of the first two targets is to constrain the reflection 
component with the best accuracy, aiming at the brightest new sources in the 
(then) latest BAT catalog. For the four targets with 20 ks exposure, the main 
goal was to obtain their X-ray spectra below 10 keV to make the BAT sample 
complete. 



and Swift J0350. 1-5019, we applied spaced-row charge in- 
jectio n (SCI) for the XIS d ata to improve the energy resolu- 
tion ((Nakaiima et al. 2008). In the spectral analysis, we also 
utilize the BAT spectra covering the 15-200 keV band, inte- 
grated for the first 9 months for Swift J0138.6-4001 and Swift 
J0601. 9-8636, and for 22 months for the rest of targets. 

3. ANALYSIS AND RESULTS 

We analyze the Suzaku data using HEAsoft version 6.3.2 
from the data products version 2.0 distributed by the Suzaku 
pipeline processing team. In extraction of the light curves 
and spectra, we set the source region as a circle around the 
detected position with a radius of 1 .5-2 arcmin, depending on 
the flux. For the XIS data, the background was taken from a 
source-free region in the field of view with an approximately 
same offset angle from the optical axis as the source. For 
the HXD/PIN data, we use the so-called "tuned" background 
model provided by the HXD team. Its systematic errors are 
estimated to be ~ 1.3% at a la co nfidence level in the 15- 
40 keV band for a 10 ks exposure dMizuno et al.ll2008l) . As 
our exposures are «20 ks or larger, we expect that the error is 
even smaller than this value. 

3.1. Light Curves 

Figures [T]and [2] show the background-subtract light curves 
of our targets obtained with the XIS and HXD in the 2-10 
keV and 15-40 keV band, respectively. To minimize any sys- 
tematic uncertainties caused by the orbital change of the satel- 
lite, we merge data taken during one orbit (96 minutes) into 
one bin. Then, to check if there is any significant time vari- 
ability during the observations, we perform a simple x 2 test 
to each light curve assuming a null hypothesis of a constant 
flux. The resultant reduced x 2 value and the degree of free- 
dom are shown in each panel. We detect no significant time 
variability on a time scale of hours for all the targets in both 
energy bands. Thus, we analyze the spectra of all the targets 
averaged over the whole observation epoch. 

3.2. BAT Spectra 

Before performing the spectral fit to the Suzaku data, we 
firstly analyze only the BAT spectra in the 15-200 keV band 
to constrain the high energy cutoff in the continuum. It 
is known that the incident photon spectra of Seyfert galax- 
ies are well approximated by a power law with an ex- 
ponential cutoff (cutoff power law model), represented as 
AE~ r exp{-E /E cul }, where A, T, E cut is the normalization at 
1 keV, photon index, c utoff energy, respectively . Here we uti- 
lize the pexrav code by Magdziarz & Zdziarski ( 1995) to take 
into account possible contribution of the Compton reflection 
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TABLE 2 
Observation Log 



SWIFT 


Start Time (UT) 


End Time 


Exposure 2 (XIS) 


Exposure (HXD/PIN) 


SCI b 


J0138.6^001 


2006 Jun()4 18:13 


Jun 05 05:00 


21.2 ks 


16.5 ks 


Off 


J0255.2-0011 


2007 Jan 23 14:54 


Jan 26 05:30 


101.6 ks 


80.6 ks 


On 


J0350.1-5019 


2006 Nov 23 02:07 


Nov 23 13:33 


19.3 ks 


15.7 ks 


On 


J0505.7-2348 


2006 Apr 01 22:12 


Apr 04 02:55 


78.6 ks 


50.3 ks 


Off 


J060 1.9-8636 


2006 Apr 13 16:24 


Apr 14 01:52 


19.8 ks 


15.7 ks 


Off 


J1628.1+5145 


2006 Apr 15 18:49 


Apr 16 09:22 


23.6 ks 


19.5 ks 


Off 



a Bas ed on the good t ime interval for XIS-0. b With/without the spaced-row charge injection for the 
XISlNakajima et al. 2008). 



Bin time: 5760 



(a) SWIFT J0138. 6-4001 



reduced t (dof) = 1.32 (6) 



10 4 2*10 4 3x10" 4x10 4 5x10" 6-1 4 

Time (s) 

Start Time 13890 19:02:29:184 Slop Time 13091 4:02:29:184 



(b) SWIFT J0255.2-001 1 



reduced (dof) = 0.61 (39) 



Time (s) 

Start Time 1 41 23 1 5:41 :34:698 Stop Time 1 41 26 5:1 1 :34:6! 



(c) SWIFT J0350.1 -501 9 



reduced % 2 (dof) = 0.39 (6) 

10 4 2x10 4 3*W 4*10 4 5x10 s 

Time (s] 

Start Time 14062 3:40:58:184 Stop Time 14062 12:40:58:184 



:n time: 5760 



(d) SWIFT J0505. 7-2348 



reduced % z (dof) = 1.46 (32) 



Time (s) 

Start Time 13827 0:02:47:184 Stop Time 13B29 3:02:47:184 



(e) SWIFT J0601. 9-8636 



reduced t (dof) = 0.21 (5) 



10 4 2x10 4 3x1 4 4x10 4 5x10" 6x1 4 

Time (s) 

Start Time 13838 17:58:31:184 Slop Time 13839 1:28:31:184 



(f) SWIFT J1628.1+5145 



+ 



-+++++++ 



reduced x z (dof) = 0.76 (9) 



10 4 2*10 4 3x10 4 4*10 4 5*10 4 

Time (s) 

Start True 13840 19:40:03:184 Stop Time 13841 9:10:03:184 



FIG. 1 . — The background subtracted light curves of XIS in the 2-10 keV band during the Suzaku observations. One bin corresponds to 96 minutes. The data 
from the XIS-0 and XIS-2 are summed for Swift J0255.2-0011 and Swift J0350. 1-5019, while those from XIS-0, XIS-2, and XIS-3 are summed for the rest. 
The numbers listed in each panel represent the value of reduced \ 2 with the degree of freedom for the constant flux hypothesis. 



component from optically thick, cold matter. The relative in- 
tensity of the reflection component to that of the incident cut- 
off power law component is defined as R = fl /2ir, where fl 
is the solid angle of the reflector (R = 1 corresponds to the 
reflection from a semi-infinite plane). 

In this stage, we assume R = or 2 as the two extreme 
cases just to evaluate the effects of the reflection. The inclina- 
tion angle is fixed at 60°. To avoid strong coupling between 
the power law slope and and cutoff energy, we fix the photon 



index at 1.9, which is a canonical value for AGNs. Table [TJ 
summarizes the fitting results of the cutoff energy for each 
target. While it is difficult to constrain its upper limit due to 
the limited band coverage of the BAT, we find the cutoff en- 
ergy must be above «100 keV in most cases. Considering 
this limitation, we fix E cut at 300 keV in all subsequent anal- 
ysis. Impacts by adopti ng a lower value of £ C u t in the spectral 
fitting are discussed in 83.51 
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TABLE 3 

Cutoff Energies (E C ut) determined by the BAT spectra 



n/2w 


SWIFT JO 1 3 8 . 6-400 1 


J0255.2-0011 


J0350. 1-5019 


J0505.7-2348 


J0601.9-8636 


J1628.1+5145 





> 317 


> 289 


> 106 


> 318 


> 116 


> 99 


X 2 /d.o.f. 


1.58/2 


4.57/6 


1.16/6 


13.58/6 


4.96/2 


2.38/6 


2 


329(>91) 


> 253 


> 95 


> 306 


> 82 


> 87 


X 2 /d.o.f. 


1.05/2 


5.21/6 


2.40/6 


15.24/6 


3.37/2 


4.61/6 


Note. — 


- The unit of £ cut is keV. 
















Bin time: 5760 






B 


"i lime: 5760 s 



(a) SWIFT J0138.6-4001 



reduced % 2 (dof) = 1.44 (6) 

10 4 2x10 4 3*10" 4*10 4 5x10 4 

Time (s) 

Start Time 13890 19:02:29:184 Stop Time 13891 5:32:29:184 



(b) SWIFT J0255. 2-001 1 



lit 



reduced (dof) - 1.07 (38) 



Time (s) 

Start Tims 14123 15:45:50:698 Stop Time 14126 5:15:50:698 



(c) SWIFT J0350. 1-5019 



reduci d % 2 (dof) = 0.86 (6) 



10 4 2-10 4 3x10" 4x10 4 5x10" 6x1 4 

Time (s) 

Start Time 14062 4:09:45:184 Stop Time 14062 13:09:45:184 



(d) SWIFT J0505. 7-2348 



tAl- 



+ tT T 



reduced* 2 (dof) = 0.69 (29) 



5x10" 10 5 1.5x10 s 2x10 5 2.5x10 s 

Time (s) 

Start Time 13827 0:23:31:184 Stop Time 13829 3:23:31:184 



Bin time: 5760 



in time: 5760 



(e) SWIFT J0601. 9-8636 



reduced % 2 (dof) = 0.39 (5) 



10 4 2x10 4 3x10" 4x10 4 5x10" 6x1 4 

Time (s) 

Start Time 13838 17:58:31:184 Stop Time 13839 1:28:31:184 



(f) SWIFT J1628. 1+5145 






+ 




+++- 




■++ +" 










reduced x 2 (dof) = 1 -04 (9) 




10 4 2x1c 


3.10 4 4x10 J 5x10 4 
Time (s) 


6* 



Stan Time 13840 19:40:03:194 Stop Time 13B41 9:10:03:134 



FIG. 2. — The background subtracted light curves of HXD/PIN in the 15^10 keV band during the Suzaku observations. One bin corresponds to 96 minutes. 
The numbers listed in each panel represent the value of reduced \ 2 with the degree of freedom for the constant flux hypothesis. 



3.3. Spectral Models 

We consider three basic models in the spectral analysis uni- 
formly for all the six targets. Our policy is to start with the 
simplest model for each target; if we find that the fit with a 
simple model does not give a physically self-consistent pic- 
ture and/or that the fit is significantly improved by introduc- 
ing additional parameters, then we adopt more complicated 
models. In all the cases, we assume a power law with an ex- 
ponential cutoff fixed at 300 keV for the incident continuum, 
as explained in the previous subsection. The Galactic absorp- 
tion, A^ al , is always included in the model (even if not ex- 



plicitly mentioned below) by assuming the hydrog en column 
density from the H I map of iKalberla et alj d2005l) . available 
with the nh program in the HEAsoft package. For absorp- 
tion, we adopt the photo-electric absorpti on cross section by 
iBalucinska-Church & McCammonl d 19921) ("bcmc") and use 
the phabs or zphabs mod el in XSPEC. Solar abundances by 
lAnders & Grevessa (119891) are assumed throughout the anal- 
ysis. 

The first (simplest) model, designated as Model A, con- 
sists of (i) a transmitted component (a cut-off power law ab- 
sorbed by cold matter), (ii) a scattered component (a cut-off 
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power law without absorption), and (iii) an iron-K emission 
line (a gaussian), represented as zphabs*zhighect*zpowerlw 
+ const*zhighect*zpowerlw + zgauss in the XSPEC termi- 
nology 6 . According to the unified scheme, in type 2 AGNs, 
the incident power law from the nucleus will be absorbed by 
a dusty torus around the SMBH, while the nuclear emission is 
partially scattered into the line-of-sight by ionized gas around 
the torus. For simplicity, we assume that the scattered compo- 
nent has the same shape of the incident power law with a frac- 
tion of / scat as first-order approximation, although, in reality, 
it often contains a number of re combination fine s produced 
by a photo-ionized plasma (e.g., ISako et aT]|2000l) unless the 
scatterer is fully ionized. As a result of reprocessing by cold 
matter in the surrounding environment (such as a torus, an 
ionized gas, and an accretion disk), an iron-K emission will be 
produced at the rest-frame 6.4 keV. Since a line from the torus 
should not resolved by the energy resolution of the XIS, we 
fix the 1 a line width of the iron-K line at the averaged value of 
the (apparent) line width of the 55 Fe calibration source at 5.9 
keV: 0.3 eV, 38 eV, 30 eV, 0.8 eV, 1.7 eV and 3.2 eV for Swift 
J0138.6-4001, J0255.2-0011, J0350.1-5019, J0505.7-2348, 
J0601.9-8636, and J1628. 1-5145, respectively, to take into 
account a possible systematic error in the energy response 1 . 

In the second model, Model B, we consider an addi- 
tional contribution of an absorbed Compton reflection com- 
ponent from optically thick matter. This emission is ex- 
pected from the inner wall of an (optically thick) torus 
and/or the accretion disk, irradiated by the incident contin- 
uum. Model B is expressed as zphabs*zhighect*zpowerlw 
+ const*zhighect*zpowerlw + zgauss + zphabs*pexrav in 
the XSPEC terminology, where the last term represents the 
reflection component (not including the direct one). In the 
pexrav model, we make the solid angle Vl of the reflector 
seen from the nucleus as a free parameter, and fix the incli- 
nation angle and the cutoff energy at 60° and 300 keV, re- 
spectively. The photon index and normalization of the power 
law are linked to those of the transmitted component. The 
absorption to the reflection component is set to be indepen- 
dent of that for the transmitted one, by considering a different 
geometry of the emission region. To check the physical va- 
lidity of Model B, we examine the equivalent width (E.W.) of 
the iron-K emission line with respect to the reflection com- 
ponent, E.W. refl , obtai ned from the fit. Theoretically, it is ex- 
pected to be 1-2 keV dMatt et al.ll 19911) . This value, however, 
tightly depends on the geometries of the torus and the accre- 
tion disk. Thus, we regard the result of Model B as valid only 
if E. W. refl = 0.5 - 2 keV but discard it otherwise. 

We finally consider the third model (Model C), where 
we assume two different absorptions with different cover- 
ing factors for the transmitted component. This reduces the 
contribution of a less absorbed reflection component com- 
pared with the case of a single absorption as assumed in 
Model B. In fact, sim ulations show tha t the absorber in the 
torus can be patchy dWada & Norman! |2002|) . resulting in 
a time-averaged spectrum that is better modeled by multi- 
ple absorptions to the transmitted continuum rather than by 
a single absorber. In the XSPEC terminology, Model C 
is expressed as zphabs*zpcfabs*zhighect*zpowerlw + 

6 the zhighect model is utilized to represent exponential cutoff at the 
source redshift. 

7 The large line widths of Swift J0255.2-001 1 and J03 50. 1-5019 are arti- 
fact caused by inaccurate calibration for the data with SCI (Matsumoto 2007; 
see http://www.astro.isas. ac.jp/suzaku/analysis/xis/ver2. 0/) 



const*zhighect*zpowerlw + zgauss + zphabs*pexrav, 

where the multiple terms of the first component (absorbed, 
partial covering model) represents the two different absorp- 
tions. 

To summarize, we can write the three models of the photon 
spectrum F (E) without the Galactic absorption as follows: 

F (E) = f c exp {-A^cr (£)}/(£) + ( 1 - / c ) exp { -iV> (£)}/ (E) 

+fsc-J (E) + G(E) + exp {-N^a (E) }C(E)+S(E), ( 1) 

where / c is the fraction of the more heavily absorbed compo- 
nent (f c = 1 in Models A and B), and are the absorp- 
tion column densities for the transmitted component (N^ = 
in Models A and B), a (E) is the cross section of photoelectric 
absorption, 1(E) = AE~ r exp{-E / E cut } is the intrinsic cutoff 
power law component, / scat is the scattered fraction, G(E) is 
a gaussian representing the iron-K emission line, C(E) is the 
Compton reflection component (C(E) = in Model A), and 
S(E) represents additional soft components (see below). 

3.4. Fitting Results 

Using the three models described in the previous subsec- 
tion, we perform spectral fit simultaneously to the spectra of 
the FI-XISs (those of 2 or 3 XISs are summed), the BI-XIS, 
and the HXD/PIN. Based on the fitting results, we chose the 
best model in the following manner. (1) First, if Models B is 
found to significantly improve the fit from Model A by per- 
forming an F-test, we tentatively adopt Model B as a better 
model, otherwise Model A. (2) Then, if Model B is found 
to be physically not self-consistent in terms of the E.W. of 
the iron-K line as explained in the previous subsection, or if 
Model C significantly improves the fit compared with Model 
A, we adopt Model C as the best model. 

After finding the best one among the three models, we 
finally include the Swift/BAT spectrum to the Suzaku data 
to constrain the photon index most tightly, assuming that 
the spectral variability is negligible between the observation 
epoch of Suzaku and that of Swift. In the fit, we fix the relative 
normalization between the FI-XISs and the PIN at 1. 1 based 
on the calibration of Crab Nebula (Ishi da et alj |2007), while 
those of the BI-XIS and the BAT against FI-XISs are treated 
as free parameters, considering the calibration uncertainty and 
time variability (in flux), respectively. The detailed results of 
spectral fit for each source is summarized below. 

Swift J01 38.6-4001 

We adopt Model B as the most appropriate model of Swift 
J0138.6^1001. We obtain (xV) = (101.14, 96) with Model 
A and (92.52, 94) with Model B, where v is the degree of 
freedom. Thus, the improvement of the fit by adding a re- 
flection component is found to be significant at >90% confi- 
dence level by an F-test, which gives an F-value of 4.38 for 
the degrees of freedom of (2, 94). Since the E.W. of the iron- 
K line with respect to the reflection component is found to 
be E.W. refl > 0.5 keV, the model is physically self-consistent. 
Here we allow the absorption to the reflect ion component N„ n 
to be a free parameter, unlike the case in lUeda et alJ 42007), 
where it is linked to that for the transmitted component. The 
basic results are essentially unchanged, however. 

Swift J 02 5 5. 2-00 11 

We adopt Model B with a apec plasma model for Swift 
J0255.2-001 1. Models A and B yield (x 2 , u) = (596.85, 304) 
and (500.75, 302), respectively, suggesting the presence of 
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FIG. 3. — The observed spectra of the 6 targets. The black, red, magenta and blue crosses represent the data of the FI-XIS, BI-XIS, HXD/PIN, and BAT, 
respectively, showing their la errors in vertical direction. The spectra of the XIS and PIN are folded with the detector response in units of counts s keV~', 
while those of the BAT is corrected for the detector area and have units of photons cirT 2 ks _1 keV~' . The best-fit models are plotted by lines, and the residuals in 
units of x are shown in the lower panels. 



a significant amount of a reflection component. We find 
E.W. refl = 0.81^ keV with Model B, which is physically 
consistent. We also confirm that applying Model C does not 
improve the fit significantly, yielding f c = 1.0 (i.e., multi- 
ple absorbers for the transmitted component are not required 
when we consider a reflection component). As positive resid- 
uals remain in the energy band below 1 keV, we add the 
apec model 8 in XSPEC, an emission model from an optically- 
thin thermal plasma with Solar abundances, over Model B. It 
greatly improve the fit, giving (x 2 ,^) = (330.38, 300). The 
temperature of the plasma is found to be kT = 0.74^} ( [q keV 
with an emission measure of n 2 V « 5 x 10 63 cm" 3 . This tem- 
perature is similar to those found in many Seyfert 2 galaxies 
irTurneretalJll997l) . 

Swift J0350.1-5019 

Model C is adopted for Swift J0350.1-5019 since we find 
E.W. refl < 0.1 keV . We obtain (xV) = (34.93, 36) with 
Model A, and (28.24, 35) with Model B. Here, we link the 
absorption to the reflection component and that to the trans- 
mitted one, as having the former as a free parameter does not 
help improving the fit. Model C gives (\ , v) = (28.24, 35) 
keV, which is now physically consistent. 



andE.W. refl = 2.8: 



http ://cxc . harvard, edu/ atomdb/ 



Swift J0505. 7-2348 

We adopt Model C with a line feature at «0.9 keV as the 
best description for Swift J0505. 7-2348. We find very little 
reflection component in the spectra of Swift J0505. 7-2348; 
Models A and B give a very similar \ 2 value of 671.3 for a 
degree of freedom of 602 and 600, respectively, with no im- 
provement of the fit by Model B. In fact, we obtain a tight 
upper limit for the reflection component of R < 0.11 (90% 
confidence level) from the fit with Model B. Thus, we apply 
Model C (here A^ efl is tied to since 7? ~ 0). This model 
yields (% 2 , v) = (636.11, 598), which is significantly better 
than Model A (or B). Furthermore, positive residuals remain 
in the 0.7-1 keV band, centered at ~ 0.9 keV. To model this 
feature, we firstly add an apec component to Model C like 
in the case of Swift J0255. 2-001 1. This is found to be un- 
successful, largely overestimating the 1-1.2 keV flux of the 
data. Rather, in some obscured AGNs, emission lines from 
Ne ions in a photo-ioniz ed plasma are obser ved around 0.9 
k eV (e.g., NGC 45 07 in IComastri etal][T998l and NGC 4151 
in lOgle et alj|2000l) . which is a likely origin of the excess fea- 
ture seen in our spectra. As these emission lines cannot be 
resolved by CCDs, we model it with a gaussian whose la 
width is fixed at 0.05 keV, referring to the result of the Seyfert 
2 galaxy NGC 4507 ( iMatt et alj|2004h . The center energy of 
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FIG. 4. — The best fit spectral model of the 6 targets in units of EFe (where E is the energy and Fe is the photon spectrum). The black, dashed red, dotted 
blue, dot-dot-dashed cyan, dot-dashed magenta curves correspond to the total, transmitted one, reflection component , scattered component, and iron-K emission 
line, respectively. The purple dashed model in Swift J0255. 2-001 1 represents the emission from an optically-thin thermal plasma, and the orange one in Swift 
J0505. 7-2348 the Ne emission lines from a photo-ionized plasma (see text). 



the gaussian is treated as a free parameter, and is found to 
be 0.87 ± 0.02keV, consistent with our assumption. We ob- 
tain (x 2 , v) = (615.08, 597), which is a significantly better fit 
compared with Model C without the line. 

Swift J0601. 9-8636 

Model B is adopted for Swift J 60 1.9 -8636. We confirm 
the results reported in lUeda et all (120071) . Models A and B 
yield (x 2 ,f) = (45.20, 37) and (32.32, 35), respectively, sug- 
gesting a significant amount of a reflection component. As 
E.W. refl is found to be 2. 10!" f keV, we find Model B to be a 
physically consistent model for this source. 

Swift J 1628.1-5 145 

Model C is adopted for Swift J1628. 1-5 145. Model A gives 
(X 2 ;^) = (574.66, 546), which is significantly improved to 
(564.23, 544) in Model B but with E.W. refl of < 0.1 keV. 
Hence, we apply Model C to the data, which yields (x 2 , v) = 
(553.88, 543) and a physically consistent value of E.W. refl = 
0.87^35 keV. 



3.5. Results Summary 

Table Q] summarizes the choice of the spectral model (A, B, 
or C) and best-fit parameters from the Suzaku + BAT simul- 
taneous fit for each target. The observed fluxes in the 2-10 
keV and 10-50 keV bands and the estimated 2-10 keV in- 
trinsic luminosity corrected for absorption are also listed. We 
find that all the six targets are heavily obscured with column 
densities larger th an > 3 x 10 23 cm" 2 . As already reported by 
lUeda et al.1 d2007h . Swift J0601 .9-8636 can be called a mildly 
"Compton thick" AGN, since it shows log Nn > 24 cm" 2 but 
a transmitted component is still seen in the hard X-ray band 
above 10 keV. We confirm that the photon indices are within 
the range of 1.6-2.0 for all the targets, and hence the fitting 
model is physically proper. 

The upper panels of Figure [3] show the observed spectra of 
the FI-XIS (black), the BI-XIS (red), and the HXD/PIN (ma- 
genta) folded with the detector response in units of counts s" 1 
keV" 1 , together with the BAT spectra (blue) corrected for the 
detector area in units of photons cm" 2 ks" 1 keV" 1 . The best-fit 
models are superposed by solid lines. The lower panels show 
the residuals in units of x (i-e., divided by the la statistical 
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error in each bin). Figure [4] shows the incident spectra with- 
out Galactic absorption in units of EFe, where the contribu- 
tion of each component in equation (1) is plotted separately; 
the black, red, blue, cyan, magenta curves correspond to the 
total, transmitted component, reflection component, scattered 
component, and iron-K emission line, respectively. For Swift 
J0255.2-0011 and Swift J0505.7-2348, the additional soft 
components are also included in purple and orange, respec- 
tively. 

For three sources Swift J0138.6-4001, Swift J0255.2- 
001 1, and Swift J060 1.9-8636, the best fit reflection strength 
parameter R(= H/2tt) > 1, although values close to unity 
are allowed within the uncertainties. Most simple geome- 
tries only allow R f« 1 to be the maximum physically plau- 
sible value. The most likely explanation for these possibly 
unphysical results is that a part of the direct emission from 
the nucleus is completely blocked by nonuniform material in 
the line of sight, reducing the observed normalization of the 
transmitted component smal ler than the true value by a factor 
of < 1/R (lUeda et al.ll2007l) . Alternatively, time variability 
can be responsible for this, if we are observing an echo of 
previously brighter phases of the AGN in the reflection com- 
ponent. Thus, we list in Table [T]the intrinsic luminosities cor- 
rected for either of these effects by multiplying R. Accord- 
ingly, the scattered fraction / scat is corrected by 1/R from the 
observed value, even though this correction does not affect 
our conclusion. 

To evaluate maximum possible errors in the spectral param- 
eters due to the uncertainty of the cutoff energy, we also per- 
form fitting by adopting E cat =100 keV instead of 300 keV for 
Swift J0350.1-5019, Swift J0601.9-8636, and Swift J1628.1- 
5145, whose 90% lower limits of £ C ut are smaller than 100 
keV as determined with the Swift /BAT spectra (§ 13. 2b . We find 
that Swift J0350. 1-5019 shows a smaller scattered fraction 
(/scat = 0.77+^1%) compared with the case of £ cut = 300 keV. 
For Swift J1628. 1-5145, we also obtain a slightly smaller 
scattered fraction (/ sca t = 0.74 ± 0.06) with stronger reflec- 
tion (R = 0.80!^) and larger absorption (N& = (87 ± 14) x 
10 22 cm" 2 ). The other parameters, including those of Swift 
J0601. 9-8636, do not change within the 90% errors. These 
results indicate that coupling of the spectral parameters with 
E cut could not always be negligible. Nevertheless, we confirm 
that it does not affect our discussion below, even considering 
the largest uncertainties as estimated here. 

4. DISCUSSION 

We have for the first time obtained broad-band spectra cov- 
ering the 0.5-50 keV band of the six "new" AGNs detected 
in the Swi'/f/BAT survey that did not have precise X-ray ob- 
servations, inclu ding the two sources already reported by 
lUeda et ail (120071) . These six targets were essentially selected 
without biases, and hence can be regarded as a representa- 
tive sample of new SwiftlBKY AGNs. The spectra of all the 
targets are uniformly described with a spectral model consist- 
ing of a heavily absorbed (log A^h > 23.5 cm" 2 ) transmitted 
components with a single or multiple absorptions, scattered 
lights, a Compton reflection component from optically-thick 
cold matter, and an iron-K emission line at 6.4 keV in the rest- 
frame, with additional soft X-ray components in two cases. 
Thanks to its good sensitivity in the 10-50 keV band and si- 
multaneous band coverage with Suzaku, we have accurately 
derived key spectral parameters including the fraction of the 
scattered component / scat , the strength of the reflection com- 
ponent R(= VL/2n), and the E.W. of the iron-K emission line. 



The first notable result is that our sample show very small 
values of / scat , < 3%, compared with a typical value of op- 
tical selected Seyfert-2 samples of 3-10% (e.g . , iTurner et al. | 
119971: ICappi et al.ll2006t iGuainazzi et~ai1l2005l) . [winter et all 
(2008bj) categorize AGNs with / scat < 3% as "hidden" AGNs, 
to which all the six sources belong. This is not unexpected, 
since the "hidden" AGNs constitutes a significant fraction 
(24%) of the uniform SwiftfBAT sample. 

Figure [5] shows the correlation between the reflection 
strength R versus scattering fraction / scat for our sample. This 
plot suggests that there are two groups even within the "hid- 
den" population defined by Winter et al. 2008b, that is, 

• /?>0.8and/ scat <0.5% 

• fl<0.8and/ scat >0.5%, 

although it is hard to conclude if this distribution is really 
distinct or continuous, given the small number of the cur - 
rent sample. Following the description bv lUeda et al.l d2007l) . 
we refer to the former group as "new type" and the other 
"classical t y pe". T he new group includes the two sources of 
lUeda et al.l (I2007T) . Swift J0138.6-4001 and Swift J0601.9- 
8636, and Swift J0255.2-001 1. The extremely small scatter- 
ing fraction indicates two possibilities (1) the opening angle 
of the torus is small and/or (2) there is very little gas around 
the torus to scatter the X-rays. Although it is difficult to firmly 
exclude the latter possibility, the global presence of a reflec- 
tion component strongly supports the former, as discussed be- 
low. 

We can constrain the geometry of the torus and viewing an- 
gles from the spectral parameters, the scattering fraction, the 
E.W. of the iron-K line, and absorption column density, based 
on a simple torus mod el. Here we utilize the calculation by 
iLevenson et al.l d2002l) . where three free parameters are intro- 
duced for a uniform-density torus whose cross-section has a 
square shape; the thickness of the torus in terms of the op- 
tical depth for Thomson scattering r, the half-opening angle 
of torus 6, and the inclination angle i. In this geometry, the 
line-of-sight column density of the torus depends on z"; when 
we see the nucleus with face-on view, the absorption becomes 
much smaller than the case of edge-on view. If the column 
density of the scatterer is constant, the opening angle can be 
connected to the scattering fraction as 



cosf? = 1 - 



/sc 



■ (1 - cos 6*o ) . 



(2) 



sc;it.() 



• New Type 
o Classical Type 



-O- 



0.2 



0.5 
f 



(%) 



FIG. 5. — The correlation between the strength of the Compton reflection 
component (R = Q/2n) and the fraction of the scattered component (/ sca i) 
for the 6 targets. Filled and open circles represent "new type" (R > 0.8 and 
/scat < 0.5%) and "classical type" (R < 0.8 and / scat > 0.5%), respectively. 
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where we normalize / scat = / scat , at 6 = 6 . Assuming / scat;0 = 
3% at 6*o = 45° as typical parameters 9 , we obtain 9 ~ 20° 
for the new type grou p. Then, comparing with Figure 2 of 
iLevenson et alj d2002l) where r = 4 is assumed, we can con- 
strain i ~ 40° for Swift J0601. 9-8636, and i ~ 20° for Swift 
J0138.6-4001 and Swift J0255.2-001 1, from the observed 
E.W. of the iron-K line (here we have neglected the contri- 
bution of an iron-K line originating from the accretion disk, 
which is often broader than the "narrow" one considered in 
the spectral model.). The absorption column densities A^h 
> iq23 .8 cm -2 an( j s t ron g reflection R > 1 are fully consistent 
with the geometry considered here. The reflection continuum 
is attributable to that from the inner wall of the tall torus and 
that from the accretion disk seen with a small inclination an- 
gle. 

Thus, it is very likely that we are seeing the new type 
AGNs in a rather face-on geometry thro ugh a geometricall y 
and optically thick torus, as discussed in lUeda et al.l d2007l) . 
An important consequence of this interpretation is that, if 
we view the same objects in a more edge-on geometry, then 
these AGNs look heavily Compton-thick (jVh ~ 10 25 cm -2 ) 
and hence the transmitted components hardly escape toward 
us even in the hard X-ray > 10 keV regime due to repeated 
scattering. Thus, similar systems to new type AGNs discov- 
ered here may still be largely missed even in the on-going 
Svw/f/BAT or INTEGRAL surveys, requiring even more sen- 
sitive observations in hard X-rays to detect them more com- 
pletely. 

For the classical type of AGNs, we obtain 9 > 30° from 
the observed scattered fraction by equation (0. Two sources 
in this group, Swift J0505.7-2348 and Swift J1628.1-5145, 
show weak E.W. of the iron-K emission line, ~ 0.05 ke V. This 
is hard t o explain by a torus wi th large r as assumed in Fig- 
ure 2 of ILevenson et al.l (120021) . which predicts E.W. > 0.1 
keV if i > 9 (i.e., the transmitted component is absorbed). 
To consider a torus with different optical depths, we apply a 
model developed by Ghis ellini et al.l (1 1994ft (see their Figure 1 
for the definition of the geometry). We find that the E.W of 
the iron-K line and absorption can be consistently explained, 
within a factor of 2, if the column density of the torus in the 
equatorial plane is ~ 10 23 " 24 cm -2 and 60° < i < 84°. The ab- 
sence of strong reflection components is also consistent with 
the small optical depth of the torus and with the "edge-on" 
accretion disk. 

To summarize, we have discovered from our first Swift- 



Suzaku AGN sample two different classes of obscured AGN 
population, called "new type" that have an extremely small 
scattering fraction (/ sca t < 0.5%) and a strong reflection com- 
ponent (R > 0.8), and "classical type" with / scat > 0.5% and 
R < 0.8. Figure [6] shows a schematic illustration of the torus 



Black Hole 




Classical Type 



New Type 



FIG. 6. — The schematic illustration of the torus geometry for the classical 
type AGNs (left panel) and the new type AGNs (right panel). 

geometry of the two types. It is likely that the new type AGNs 
are deeply buried in an optically thick torus (A^h ~ 10 25 cm" 2 ) 
with a small opening angle 9 ~ 20°, and are viewed in a face- 
on geometry. By contrast, the classical type AGNs have an 
optically thin torus (A^h ~ 10 23-24 cm" 2 ) with a larger open- 
ing angle 9 > 30° viewed in a more edge-on geometry than 
the new type. A significant number of new type AGNs with 
an edge-on view from us may still be missing in the current 
all-sky hard X-ray surveys. The presence of two classes of 
obscured AGNs implies that the understanding of AGN is not 
as simple as assumed in the population synthesis models of 
the cosmic X-ray background. Given the small number of our 
sample, however, it is not clear whether there are two "dis- 
tinct" types of obscured AGNs or the distribution is contin- 
uous. More systematic follow-up observations with Suzaku 
of new Svw/f/BAT AGNs are important to establish the true 
properties of local AGNs. 
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9 The value 0q = 45° corresponds to the number ratio of obscured AGNs 
to unobscured ones of about 3, which is consistent with the observations of 



local AGNs. 
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Best-fit Spectral Parameters 
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1.0 x 10~ 12 


9.0 x 10" 12 


(12) 


F10-50 (ergscnr 2 s~') 


3.4 x l(T n 


4.3 X lO" 11 


5.5 x 10~ 12 


2.9 x lO" 11 


1.1 x 10-" 


3.6 x 10-" 


(13) 


L2-10 (ergss -1 ) 


2.4 x 10 43 


5.1 x 10 43 


1.4 x 10 43 


1.5 x 10 44 


8.9 x 10 41 


1.5 x 10 44 




X 2 /d.o.f. 


94.42/97 


335.97/307 


27.28/39 


622.38/604 


35.48/38 


558.76/550 



r 
< 



NOTE. — (1) The hydrogen column density of Galactic absorption by Kalberla et al. 1 2005). (2) The line-of-sight hydrogen column density for the transmitted component. For the double covering model (Model 

C), that for more heavily absorbed component is given. (3) The smaller line-of-sight hydrogen column density for the transmitted component in Model C. (4) The normalization fraction of the more absorbed 

component to the total transmitted one in Model C. (5) The power-law photon index. (6) The fraction of the scattered component relative to the intrinsic power law,corrected for the transmission efficiency of 1/i? ^ 

if R > 1 (see text). (7) The center energy of the iron-K emission line at the rest frame of the source redshift. (8) The observed equivalent width of the iron-K line with respect to the whole continuum. (9) The o 

line-of-sight hydrogen column density for the reflection component. (10) The relative strength of the reflection component to the transmitted one, defined as R = Q/2n, where Q is the solid angle of the reflector 

viewed from the nucleus. (11) The observed flux in the 2-10 keV band. (12) The observed flux in the 10-50 keV band. (13) The 2-10 keV intrinsic luminosity corrected for the absorption and the transmission ^ 

efficiency of 1/R if R > 1. The errors are 90% confidence limits for a single parameter. ^ 

a An additional emission from an optically-thin thermal plasma with Solar abundances is required, modelled by the apec code with a temperature of kT = 0.74+9, "„ keV and an emission measure of 5 X 10 63 cirr 3 tt) 

b " 
(see text). An emission line feature at ~ 0.9 keV is required, probably that from Ne ions from a photo-ionized plasma (see text). H 
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n 



Co 



